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ABSTRACT 

We investigate, through 2D MHD numerical simulations, the interaction of a 
uniform magnetic field oblique to a moving interstellar cloud. In particular we explore 
the transformation of cloud kinetic energy into magnetic energy as a result of field 
line stretching. Some previous simulations have emphasized the possible dynamical 
importance of a "magnetic shield" formed around clouds when the magnetic field is 
perpendicular to the cloud motion ( [Jones et al. 1996 , Miniati et al. 199~£| ). It was not 



clear, however, how dependent those findings were to the assumed field configuration 
and cloud properties. To expand our understanding of this effect, we examine several 
new cases by varing the magnetic field orientation angle with respect to the cloud 
motion (#), the cloud-background density contrast, and the cloud Mach number. 

We show that in 2D and with 6 large enough, the magnetic field tension can become 
dominant in the dynamics of the motion of high density contrast, low Mach number 
clouds. In such cases a significant fraction of cloud kinetic energy can be transformed 
into magnetic energy with the magnetic pressure at the cloud nose exceeding the ram 
pressure of the impinging flow. We derive a characteristic timescale, r ma , for this 
process of energy "conversion" . We find also that unless the cloud motion is highly 
aligned to the magneitc field, reconnection through tearing mode instabilities in the 
cloud wake limit the formation of a strong flux rope feature following the cloud. Finally 
we attempt to interpret some observational properties of the magnetic field in view of 
our results. 
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1. Introduction 



Magnetic fields represent an important component of the interstellar medium (ISM) and 
in many cases cannot be neglected when studying the processes taking place there. Numerous 
observations have been carried out and established important properties of the galactic magnetic 
field (e.g., [Vershuur 197C , Heiles 197C , Zeldovich, Ruzmaikin fc Sokoloff 1983 ). Faraday rotation 
measures ( pardner Whiteoak 1965 , Morris & Berge 1964, Valee Kronberg 1975| ) as well 
as polarization measures ( pVlathewson fc Ford 197"c| ) have established the existence of a large 
scale component of the galactic magnetic field, whereas optical polarization of starlight and 
polarization of the background synchrotron radio continuum have been used to determine its 
geometry ( Spoelstra 1977| , Gomez De Castro, Pudritz fc Bastien 1997] ) . Measurements of the 
Zeeman splitting of the hyperfine 21 cm line of neutral hydrogen (H I) have revealed a magnetic 
field strength of several /iG for many interstellar clouds ( |Vershuur 1970] 1989, [Heiles 1989| , Myersj 
et al. 1995). In addition the magnetic field is often observed to be dynamically relevant, with 
its energy being comparable to the cloud kinetic energy and, for self gravitating clouds, to their 



gravitational energy, too ( Heiles 198E , Verschuur 198£ , Myers et al. 1995 ). Since its importance 
was recognized the magnetic field has been investigated in more and more detail. The recent 
development of accurate magnetohydrodynamical (MHD) codes has provided an important new 
tool to compare and test analytical works. MHD numerical codes also can investigate beyond the 
necessary simplifications required to reduce an otherwise overwhelming mathematic complexity of 
the problems under study. 

Magnetic fields may play an especially important role in the dynamics of cloud motion 
through the ISM. Mac Low et al. (1994), for example, carried out 2D numerical simulations of the 
interaction between a cloud a few times denser than the background medium and a plane shock 
wave and showed that even a modest magnetic field aligned with the direction of the shock motion 
can help to stabilize the cloud against disruptive instabilities. This results from stretching of the 
magnetic field into a wake that forms behind such a cloud from its motion through the postshock 
plasma. That significantly reduces the vorticity shed by the cloud, which helps to stabilize it. 
Jones et al. (1996) demonstrated the same effect for clouds set into supersonic motion through 
the ISM. The latter paper also examined in detail the magnetic field behavior during supersonic 
cloud motion perpendicular to the prevailing magnetic field. In that instance Jones et al. (1996) 
found that the magnetic field plays an even more dramatic role, because field lines caught on 
the "nose" of the cloud become highly stretched and form a "magnetic shield" with a magnetic 
pressure at least comparable to the gas ram pressure. That also effectively quenches disruptive 
Rayleigh- Taylor and Kelvin- Helmholtz instabilities (KHIs), so that the clouds maintain their full 
integrity. In addition, this enhanced magnetic pressure at the front of the cloud hastens the cloud's 
deceleration with respect to the ISM. They found this result even when the initial magnetic field 
is quite modest (i.e., wherever (3 = p g /pb < 100, for a Mach 10 cloud). 

Most recently, Miniati et al. (1997, 1998) have used 2D simulations to study head-on collisions 
between mildly supersonic clouds in the ISM. Miniati et al. (1998) found a dramatic influence of 
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the magnetic field when it is perpendicular to the motion of the clouds. If the magnetic shield 
seen by Jones et al. (1996) forms before the collision takes place, it is seen to act as an almost 
elastic bumper between the two clouds, preventing their direct contact and also their disruption 
during the interaction. 

Those results clearly indicate a need to understand more completely the interactions between 
cloud motion and magnetic field. There are several simplifications in the previously mentioned 
studies that need to be relaxed along the way to that understanding. The two most obvious are 
geometrical. All of the above studies were 2D in the sense that they assumed symmetry in one 
direction perpendicular to the cloud motion. The second is the assumption that the magnetic field 
was either exactly aligned with cloud motion or perpendicular to it. 

We are tackling the fully 3D study of this problem in a different paper in preparation (Gregori 
et al. 1998). In the transition from 2D to 3D, several new physical aspects appear, adding 
enormous complexity to the investigation process. More realistic 3D MHD numerical simulations, 
although nowadays accessible through large computers, still remain much more expensive than 2D 
ones. Extending a computational box of size N x x N y to three dimensions with size N x x N y x N z 
implies an increase in the computational cost (execution time) of the calculation by a factor about 
3 /2 N z . Since generally N x>yjZ should be at least several hundred to properly capture the dynamics, 
the tradeoff usually needed to carry out 3D simulations is a reduction of either the computational 
box size or the resolution. There is a lower limit on the box size set by the influence of the 
boundary conditions. The latter are particularly important in this case because, as we shall see, 
the evolution of the magnetic energy is significantly altered by the outflow of regions of enhanced 
magnetic field. If, on the other hand, the resolution is reduced by a factor r > 1, then the cost of 
a 3D simulation with respect to a 2D one is, instead, 3 /2 N z /r A . For example a 3D calculation 
carried out on a 256 3 grid "costs" 48 times the same 2D calculation on a 512 2 grid, compared to 
an increased cost factor of 768 on a 512 3 grid. Memory and data storage factors are comparable. 
Even the smaller factor represents a major increased investment. In addition visualization and 
3D analysis of data are much more difficult. For these reasons 2D calculations are still very 
worthwhile, especially when an exploratory parametrical study is to be carried out. This is the 
case in the present circumstance, where we are investigating the role, in the cloud-magnetic field 
interaction, of several factors such as the Mach number (M), the angle between the cloud motion 
and the magnetic field lines (9) and the cloud-intercloud density contrast (%). 

Since the magnetic bumper seen in previous studies involves perpendicular fields with clearly 
dominant features, the objective of this paper is to understand how nearly aligned the magnetic 
field can be with the cloud motion and still form an effective magnetic bumper. That leads us 
naturally also to explore more generally the exchange between kinetic and magnetic energies 
during the motion of a 2D cloud through a conducting background. It is important for this 
particular purpose to question the validity of the 2D approximation. 

First, in a 2D cartesian geometry we are limited to represent cylindrical clouds with third 
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axes perpendicular to the computational plane. Since ISM clouds are observed to take many 
forms including filaments and sheets, in this respect it is still appropriate to investigate the 
simple 2D case of a transverse cylinder. However, in 2D, instabilities along the symmetry axis 
are inhibited. This means that our calculations are reliable only up to the characteristic time for 
those instabilities to occur on a spatial scale comparable to the typical cloud size. Since the most 
relevant instabilities for dense clouds in the ISM develop on a timescale r cr = 2 x 1 / 2 /M t cs (t cs 
is the sound crossing time of the cloud, see for details), we will restrict the investigation of the 
evolution of the various cases to t ^ several r cr . 

In addition it is often remarked that in a realistic 3D situation the magnetic field lines 
would be free to slip by the cloud sides, therefore escaping from the line stretching process. 
What happens in reality is quite complex and involves not only the flow of the magnetized 
gas around the cloud, but also the backreaction of the cloud to that flow and the embedded 
field. In Miniati et al. (1998) we addressed this issue briefly, pointing out the importance of the 
deformations undergone by the cloud during its motion through a magnetized gas, and their role 
in the development of the magnetic bumper. This issue will be investigated more thoroughly in 
our forthcoming paper (Gregori et al. 1998). There we shall also address some specific dynamical 
and morphological properties related to the development of 3D instabilities. In particular, 3D 
Rayleigh- Taylor instabilities (RTIs) and the consequent escape of magnetic flux tubes from the 
stretching region. In fact, the full development of RTI in the direction perpendicular to both the 
magnetic field and the cloud motion, suppressed in our 2D simulations, is able to break up the 
cloud after several r cr , thus freeing the magnetic flux tubes from the stretching mechanism. This 
process was discussed in the context of Planetary Nebulae by Soker and Dgani (1997), Dgani and 
Soker (1998) and Dgani (1998). In the meantime, in § |4.5| we will present some preliminary results 
about the growth of the magnetic energy in an analogous 3D simulation of a cylinder moving 
through a transverse magnetic field. In support of the present investigation we mention that those 
3D results are in qualitative agreement with our 2D numerical study. 

The paper outline is as follows: §^ provides a brief background review relevant to this paper. 
In §H| the numerical setup is described. In §|] we present the results of our calculations and 
comment on their physical and astrophysical implications. In §[5] we attempt to interpret some 
recent observational results based on our findings, which are summarized in 



2. Hydrodynamics and Magnetohydrodynamics of Single Cloud Motion 



A cloud set into motion within a lower density medium produces new structures in both 
the cloud and the background medium. Several features develop, each one characterized by a 
distinct timescale. In this section we will briefly review those features most relevant to the current 
study, in order to facilitate the understanding of the issues raised in this paper (for further 
details see e.g., Murray et al. 1993; Jones et al. 1994, 1996; Schiano, Christiansen & Knerr 1995; 
Malagoli, Bodo & Rosner 1996; Vietri, Ferrara fc Miniati 1997| ). The speed of the cloud (v c ) is 



- 5 - 



usually expressed in terms of the Mach number M = v c /c s i, where c s i is the sound speed in the 
unperturbed background medium. For standard H I clouds, usually M > 1 (Spitzer 1978); that 
is, the motion is supersonic. In purely gasdynamic bow shock forms in front of the cloud 

and reaches a constant standoff distance from it on the timescale Tb s ~ 2R c /v c = 2r c where R c is 
the cloud radius and r c is the cloud crossing time. This increases the gas density on the nose of 
the cloud and leads to a pressure there comparable to the ram pressure piv 2 . In response to this, 
a "crushing" shock forms and propagates through the cloud with relative speed v cs ~ v c /x 1 ^ 2 , 
where the density contrast parameter, x = Pc/Pi, is the ratio of the cloud and intercloud medium 
densities. The timescale over which this shock crosses the cloud is 

v cs v c M 

where t cs = R c /c s i is the cloud sound crossing time. After a time ~ r cr the cloud structure 
has been significantly modified and is far from uniform. Indeed, as the cloud plows through the 
background gas, KHIs and RTIs develop on its boundary layer. The slowest growing but most 
effective modes in disrupting the cloud are those characterized by the largest wavelength A ~ R, 
and for j(> 1 their typical growth time is roughly the same as r cr . Finally, as the ram pressure 
carries out work on the cloud, the cloud velocity relative to the background medium is reduced. 



This takes place on the timescale roughly given by (Klein, McKee & Colella 1994) 



3 X R Rx x r (00 , 
Av c Mc si M K 1 

The propagation of a cloud through a magnetized medium adds new parameters to the 
characterization of the cloud evolution. One is the strength of the magnetic field, expressed in 
terms of 

p B 7 V M J V ' 

where p g and pb = B 2 /8ir are the gas and magnetic pressure respectively and Ma = v c / (B/y/4wp) 
is the Alfvenic Mach number. The other is the field orientation, which in 2D calculations with 
B z = is determined by the angle 9 between the initial cloud velocity and the magnetic field lines. 

For a large /3, the beginning of the MHD cloud evolution is similar to the HD case. Over 
time, new features develop in response to the magnetic field. As mentioned earlier, numerical 
investigations have shown significant differences between the cases of cloud motion parallel 
{aligned) or perpendicular {transverse) to the initial uniform magnetic field ( |Jones et al. 1996| , 
[Mac Low et al. 1994 ). MacLow et al. (1994) studied the MHD interaction of a high Mach number 



(M=10) planar shock wave with a non radiative low density contrast {x = 10) spherical interstellar 
cloud. They investigated both aligned and, for one cylindrical cloud case, perpendicular field 
geometries, although they mostly concentrate on the former geometry. On the other hand, Jones 
et al. (1996) studied the propagation of an individual supersonic cylindrical cloud characterized 
by the same parameters {M = 10, x = 10), through a magnetized medium. In both sets of 
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computations the initial magnetic field was uniform throughout the grid. The studies found 
in common agreement that, in the aligned case, the field lines threading the cloud are pulled, 
stretched and folded around the cloud forefront, eventually reconnecting and generating new 
and more stable flux tubes in this region. Here the magnetic field never becomes dynamically 
dominant, although it plays a relevant role in stabilizing the flow around the cloud. On the other 
hand, low pressure in the cloud wake and stretching of the field lines by the wake flow cause the 
magnetic field to prevail {(3 < 10) in this region ("flux rope") behind the cloud ( [Mac Low et 
al. 1994 H 

ones et al. 1996|) . In the shocked cloud problem ( pVlac Low et al. 1994| ) the magnetic 



field in the flux rope is amplified up to the postshock thermal pressure, whereas in the isolated 
cloud calculations ( pones et al. 1996| ) such amplification is smaller. 

When the cloud motion is transverse to the magnetic field, field lines start to drape around 
the cloud surface, being compressed and, more importantly, stretched around the cloud nose. 
In particular, because of stretching of the field lines, the magnetic field in this region becomes 
eventually very strong, generating a stabilizing magnetic shield region of small values of (3 flJonesj 
et al. 1996| ), The presence of the magnetic shield also significantly alters the outcome of cloud 



collisions in the perpendicular field case (Miniati et al. 1998). Indeed, if the magnetic shield is 
strong enough, it acts like an elastic "bumper" reversing the clouds motion before they contact. 

Both of the studies involving perpendicular fields flJones et al. 1996| , [Mac Low et al. 1994j) 
found that the amplified field on the cloud nose becomes comparable to the ram pressure of the 
impinging flow, which is of the same order as the thermal pressure of the postshock gas. We shall 
show below that those conclusions were implied by their choice of the cloud parameter {\ = 10 
and M = 10). Then the question raised is how one may more generally relate the magnetic field 
properties to the initial conditions. Such question is quite complex and will be addressed in some 
more detail in §4.7. We shall emphasize though that the enhancement of the magnetic energy is 
in general due primarily to the work done on field lines through stretching by the moving cloud. 
Therefore, in principle, the cloud initial kinetic energy, 1 /2 M c v\, is the source and the upper 
limit for the magnetic energy. Nevertheless, as we will see in the following, a cloud is not always 
in the right setting for its kinetic energy to be transformed into magnetic form. However, when 
such conditions exist, the field enhancement maybe larger than previously estimated by a factor 
in principle up to ^ x which, for high density clouds, makes a significant difference. 



3. Numerical Setup 

To explore these issues we have carried out numerical calculations with an ideal MHD code, 
based on the conservative, explicit, second-order TVD method described in detail in Ryu & Jones 
(1995), Ryu, Jones, & Frank (1995) and Ryu et al. (1998). Recent measurements have shown that 
H I diffuse clouds are often characterized by high values of the electron fraction x e ~ 10 -2 ' 7 — 10 -4 ' 9 
(Myers & Khersonsky 1995). Therefore, both the kinetic Reynolds number, Re {=vrjv\ and the 
magnetic Reynolds number, ReM {=vI/um)i are expected to be large for velocities v ~ a few km 
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sec 1 and spatial scales t ~ 1 pc (Myers & Khersonsky 1995). The resulting ambipolar diffusion 
timescale is given by 



u/ " \10cm 3 J \pcj \pG J \10 -3 

meaning that the magnetic field should not decay through ion-neutral drift during the cloud 
evolutionary time considered here. Ambipolar diffusion can also affect the shock wave structure 
(Draine &: McKee 1993, and references therein). This is particularly important inside neutral 
clouds where the ionization fraction x e can be as small as ~ 1CP 5 . As a result the "shock 
transition", across which both ions and neutrals undergo the ordinary "jump conditions", can 
become so long that radiative cooling is relevant there (e.g., Draine & McKee 1993). However, 
this process only affects the thickness and structure of the shock wave, but not the overall jump 
conditions, which are set by conservation equations. These jump conditions are always correctly 
calculated by any proper TVD codes. Other processes like molecular thermal dissociation and 
radiation emission which might occur in such modified, radiative shocks are not in the scope of 
the present study. It is, therefore, appropriate to employ an ideal MHD code for investigating 
the problems of this paper. In such codes electrical resistivity is mimicked through numerical 
approximations, allowing such important real MHD effects as reconnection. 

We have used the multidimensional, Cartesian version of the code (Ryu, Jones, & Frank 
1995) with a new "Method of Characteristics, Constrained Transport" scheme for preserving 
the V • B = condition as described in Ryu et al. (1998). The computational domain is on 
the xy plane and both magnetic field and velocity z-components are set to zero. The resolution 
characterizing all the cases listed in Table |] is 50 zones per cloud radius. The dimensions of the 
grid are given in the fifth column of the same Table. When exact mirror symmetry exists across 
the X-axis (Cases 1, 3, 4, 8 and 9), only the plane y > has been included in the computational 
box, thus halving the grid size. Inflow conditions are set on the left boundary, whereas top and 
right boundaries are always open. The bottom boundary is also open except for Cases 1, 3, 4, 8 
and 9 for which it is reflective. 

Except that in one case (Case 9), in which radiative losses are being taken into account, we 
generally assume an adiabatic flow with index 7 = 5/3 (p = [7 — l]e). The clouds are always 
initially uniform and in pressure equilibrium with the background medium. The cloud density 
p c is greater than the background pi by a factor \ = Pel Pi- We set our units so that the initial 
pressure po = 3/5 and pi = 1, giving the background sound speed c s i = (jpo/p) 1 ^ 2 = 1- The 
initial magnetic field is also uniform and its strength is conveniently expressed by the parameter 



(3 (see eq. |2-3|) . We assume (3 = 4 in all the calculations. Since the thermal pressure (po) has 
not been specified in any physical units, the assumed value for (3 does not imply any particular 
strength of the magnetic field. However, for typical conditions in the ISM for which po/ks = 1600 
K cm -3 , (3 = 4 corresponds to B ~ 1.2/zG. Self-gravity has been neglected throughout our 
calculations. We have assumed very simple initial conditions for our simulations. In particular, 
several quantities, such as pressure and magnetic field, are initially set as uniform, which is 
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certainly not the case for the ISM. However, these approximations are not too important for the 
objective of our study. In particular the supersonic motion of the clouds is not affected by the 
initial pressure balance. Moreover, the dynamical interaction of clouds and magnetic field (with 
the deceleration of the former and the local amplification of the latter) can only be significant 
for the large scale component of the magnetic field. In fact magnetic field lines can be stretched 
only if their extent is much larger than the cloud size. On the other hand, local inhomogeneities 
would be quickly advected past the cloud causing only minor deviations. Thus, as long as there 
is an underlying "uniform" component to the ISM magnetic field on the scale of a few parsecs, 
the presence of local irregularities will not substantially alter the behavior expected. Interstellar 
polarization observations (e.g., Mathewson Ford 1970| ) and Faraday rotation measures (Gardner 



fc Whiteoak 1963| , Morris & Berge 1964, Valee fe Kronberg 1975| ) show clearly the existence of 



such an underlying uniform field. 

With Case 9 in Table p] we intended to investigate the possible effect of radiative losses on the 
cloud interaction with the magnetic field. The radiative cooling function we have used includes 
free-free emission, recombination lines, and collisional excitation lines as well as heating terms 
provided by collisional ionization and ionization by cosmic rays. Its properties are fully described 
in Ferrara & Field (1994) whereas its implementation is detailed in Miniati et al. (1998). Since 
the cooling function is not crucial for the issues related in this paper, we refrain from going into 
further details and refer to the above references for a more complete description. 

The clouds are initially set in motion with a Mach number M = v c /c s i, where c s j is the 
sound speed in the inter-cloud medium. As shown in Table |l] different Mach numbers have been 
employed in this work. In addition, we explore several values of 9 = arccos[v c • ~B/(v c B)], ranging 
from the aligned case (9 = 0°) to the transverse case (9 = 90°), as well as several values of x, as 
listed in Table [j]. 



4. Results 

4.1. Line Draping in a Oblique Magnetic Field 

It was clear from previous numerical studies mentioned in §|2| that the orientation of the 
magnetic field lines with respect to the cloud motion plays a key role in the interaction between 
the cloud and the magnetic field. It is important, therefore, to investigate in more detail the 
general case of cloud motion making an angle 9 with the initial uniform field (see Figure [l]). 
Intuitively we expect that the transverse field scenario should be the most common one. In fact, 
for any 9 > 0, as the cloud moves through the magnetized background and pushes the field lines 
at its nose, the field lines would tend to wrap around the cloud. We have tested this expectation 
by performing numerical simulations involving various choices of 9. Some results are presented 
in Figure |2|, which shows, for the two cases 9 = 45° and 10° the density distribution with the 
magnetic field lines superimposed. Each time corresponds to t = 13t cs = 19.5r c and 15r cs = 22.5r c 
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(t c = R c /v c = t cs /M) respectively. We emphasize that if the field lines are frozen in the plasma, 
then it takes about 2r c for them to be transported across the cloud diameter by the background 
flow and to start being stretched. 

The first thing we note is that in both cases {9 = 10° and 45°) the magnetic field lines drape 
around the cloud, similarly to the transverse field case ( [Jones et al. 1996 ). However it is clear 



that the magnetic field is significantly stronger for 6 = 45° than for 9 = 10°. Even though the 
field lines drape around the cloud, in order for cloud motion to substantially amplify the magnetic 
field and, therefore, generate the magnetic shield, a significant amount of line stretching needs 
to occur. This need not take place for all 9 > 0. Therefore, in the following we will attempt to 
identify more clearly the range of angles 9 over which stretching of the field lines turns out to be 
efficient. In addition, other factors are important. In particular the cloud could be decelerated 
by the ram pressure of the background medium or disrupted by instabilities before the magnetic 
field lines are substantially stretched. These processes are regulated by the cloud Mach number 
(M) and density contrast( %). It is, therefore, also in the objective of this study to determine the 
conditions, in terms of M and %■> that allow the formation of the magnetic shield. 



4.2. Magnetic & Kinetic Energy Evolution 

In order to explore these issues we performed a variety of numerical simulations, characterized 
by the different parameters listed in Table [l|. These simulations include: Cases 1-3 (set 1), 
corresponding to low density (x = 10) but high Mach number (M = 10) clouds, with 9 = 0°, 45°, 
and 90° respectively; Cases 4-9 (set 2), represent the opposite pairings of high density (x = 100), 
but low Mach number (M = 1.5) clouds, with 9 = 0°, 10°, 30°, 45°, 90° and 90° respectively. Case 
9 is identical to Case 8 except for being radiative. Finally Cases 10 and 11 (set 3) include clouds 
with x an d M values that are both either small (x = 10, M = 1.5) or large (x = 100, M = 10) 
with 9 = 45°. Figure || shows the time evolution of the magnetic energy increment normalized to 
the initial cloud kinetic energy (left panels) and of the relative kinetic energy decrement (right 
panels) for cases of set 1 (top), 2 (middle) and 3 (bottom) respectively (see caption for details). 
Time is espressed in units of the sound crossing time, which only depends on the cloud size and 
the background sound speed, but not on the parameters that differenciate the various cases listed 
in Table |l]. It is apparent that in most of the cases the magnetic energy increment amounts to 
an important fraction of the cloud kinetic energy. This is especially true for set 2, which shows 
that even for 9 = 30° that fraction is around 15%. On the other hand, the right hand plots show 
that the evolution of the kinetic energy is affected by the interaction with the magnetic field. For 
all cases in set 1 and 2 clouds are respectively characterized by the same parameters: it appears 
that the stronger the cloud interaction with the magnetic field (measured as the increment in the 
magnetic energy), the quicker the kinetic energy decrement. In particular for Cases 6, 7, 8 and 
9 of set 2 the differences in the kinetic energy roughly correspond to differences in the magnetic 
energy. It is clear also that the presence of radiative losses is no relevant the dynamical interaction 
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of the cloud with the magnetic field. We point out (and we shall justify it later) that flattening of 
the magnetic energy growth is in general due to both cloud deceleration and outflow of enhanced 
field regions. 



4.3. Magnetic Field Reconnection in the Cloud Wake 

In the following subsections we shall investigate and calculate in more detail the growth 
rate for the magnetic field strength (§4.4 and §4.6). Before that, we intend to explore another 



important issue concerning the conditions of the magnetic field in the cloud wake. Recall that for 
9 = a strong magnetic rope developed there in the previous simulations. We see that even at 
very small angles 9, because the field lines tend to drape around the cloud shape, the topology 
of the magnetic field in the cloud wake is basically that of the transverse field case. As a result, 
the post-cloud "flux rope" only occurs in very special cases and is mostly the outcome of the 
artificially exact alignement between the cloud motion and the magnetic field lines. In all the 
other cases we have carried out, magnetic field lines, converging in the cloud wake from the 
above and below, form an antiparallel region separated by a thin current sheet in the wake. This 
configuration is classically unstable to resistive, "tearing mode" instability. In fact, the current 
sheet in the transition layer between the field lines of opposite direction eventually becomes much 
thinner than its width (Biskamp 1993, p. 152 for details). The occurrence of the instability breaks 
up the transition current sheet into line currents and several closed field line loops arise as field 
lines reconnect across the sheet (e.g., Melrose 1986, p. 151). This process is clearly illustrated in 
Figure [|, which shows a sequence of three panels representing the evolution of the magnetic field 
at t/r cs = t/(Mr c ) = 12, 13 and 14 respectively, for Case 5 (6 = 10°, M = 1.5, x = 100)- In th e 
top panel the field lines in the transition layer appear distorted. Such distortions lead later on to 
the formation of field line loops there, as shown in the next panel corresponding to t = 13r c . These 
loops have formed through reconnection processes and they provide the signature of the tearing 
mode instability (Melrose 1986, p. 151). The overall process responsible for the modification of 
the magnetic field topology is commonly referred to as "tearing-mode reconnection" (e.g., Melrose 
1986). At t = 14t c small field line loops still exist, although they are annihilating and are finally 
being convected away from the cloud. 

It usually appears that when the magnetic field is oblique with respect to the cloud motion, 
more tearing-mode reconnection activity takes place in the cloud's wake, as compared to the 
exactly transverse field case. We attribute this property to the broken symmetry across the X-axis, 
which makes the flow in the wake more irregular, thus stimulating the onset of the instability. 
Moreover, the uneven magnetic tension caused by the broken symmetry of the magnetic field 
perturbs further the motion of the cloud. In fact it induces a vortical motion of the cloud body, 
which affects annihilation processes of the magnetic field inside the cloud itself. In addition, 
since the magnetic shield develops asymmetrically on the cloud nose, the cloud tends to drift 
perpendicular to its initial motion, toward the direction of lower magnetic pressure. 
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4.4. Analysis of the Simulations 

It is common to compare the magnetic energy in a cloud to its kinetic and/or thermal energies 
(Heiles 1989, Vershuur 1989, Myers & Khersonsky 1995, Myers et al. 1995), in order to find 
insightful trends among the various physical quantities of the ISM. Moreover, if the generation 
of local irregularities in the interstellar magnetic field is in part related to the motion of dense 
clouds (see §||), and in particular is generated at the expense of their kinetic energy, then such 
comparison may turn out to be relevant for the overall energy budget of galaxies. 

The equation describing the evolution of the magnetic energy Eb in a volume V due to field 
line stretching is ( phandrasekhar 1961 ) 



= T- I B ^ dV ~ 2£**Mb = 2^ sinO, (4-1) 
ot 4tt Jv oxj R c t c 

where sin# takes into account the field orientation. This generalizes eq. 10 in Jones et al. (1996). 
If Eb and Ebo are the magnetic energy in the region where the stretching of the field lines takes 
place at an arbitrary time t and at t = respectively, then an approximate solution to equation 
Tlis: 

E B = E B q e 2tsiae ' Tc . (4-2) 



Equation |4-2| describes an exponential growth of the magnetic energy with a timescale as short as 
T c /sin#. This approximation holds true only for the very beginning of the evolution. Depending 
on initial field strength, nonlinear effects {i.e., from flow modifications) soon become important 
so that a different approach must be employed. Starting from eq. |4-1| , we assume that the 
magnetic energy growth is initially exponential. A simple way to modify this for nonlinear growth 
is to begin from the representation of the exponential function as a series of power laws, each 
with its appropriate coefficient. We may express the saturation of field growth by saying that 
the back reaction of the magnetic field suppresses the contribution of high order terms that 
dominate the exponential function for large arguments. On the other hand, our simulation data 
show that the relative magnetic energy enhancement (AEb/Ebo) is faster than implied by a 
first order approximation (linear growth with time), but is well described by a power law with 
index 2 > m > 1 (see the next paragraph). The physical meaning of a power law behavior can, 
at least in part, be inferred from eq. |4-1| . Specifically it implies that the velocity shear dv x /dy, 



which in eq. 4-1 is represented as a constant ~ v c /R c , actually decreases with time like ~ t^ 1 
once the magnetic field begins to modify the flow. That behavior makes intuitive sense if one 
considers that the magnetic field affects the direction of the gas flow around the cloud before the 
cloud velocity itself is decreased. Thus, the effect is probably mostly due to the growing size of 
the magnetic shield, the region where field line stretching occurs, and an associated broadening 



of the shear layer around the cloud. That is, the scale represented by R c in eq. 4-1 apparently 



increases linearly with time. We derive an estimate of the magnetic energy growth by expanding 



the exponential function in eq. |4-2| and retaining a single (presumed dominant) term of the power 
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law series, with index m to be determined: 

Ebo E B o ^ ' V t cs J mV 

where in the last expression we used t c = r cs /M. 

Eq. was tested for all cases listed in Table ||. We plotted \og(AEB / Ebo) versus log(i/r cs ). 
If these two quantities are related as described by eq. (O, then the plots should resemble 



/ Eb \ 2 m t 

log I 1 I ~ log — - + m log M + m log(sin 9) + m log — 

\E B o J ml t cs 

= A(m,M,9) + m log— (4-4) 

1~cs 

and their slope should correspond to the index m. Log plots of the magnetic energy increment 
(this time normalized to the initial magnetic energy) as a function of time are shown in Figure ||. 
Left plots are log analogous of those in Figure || (left panels) and refer to the magnetic energy 
over the whole grid. On the other hand, right plots correspond to the magnetic energy in a region 
of size R c x 2R C and placed before the cloud's nose. As we can see, the growth of the magnetic 
field is much stronger in the latter case than in the former, so that the magnetic energy density is 
much higher at the cloud nose than anywhere else. However, because that is only a small region, 
its contribution to the total increase in the magnetic energy is not dominant. In other words, most 
of the exchange of kinetic and magnetic energy takes place on a more extended region. 

In general, in both cases the power law description of the magnetic energy growth with 
time applies only for 9 > 30°. So, we will not consider cases with 9 < 30° any further. Strictly 
speaking, only cases belonging to set 2 (and with 9 > 30°) follow the power law behavior. Curves 
corresponding to all other cases exhibit a turn over at some point in their evolution. That is 
particularly evident in set 1. There are two main reasons for that behavior and both are related 
to the fact that cases of sets 1 and 3 are characterized by high Mach numbers (M) and/or low 
contrast density (x)- In fact such parameters determine a quick deceleration of the cloud by the 
background ram pressure (§|2|). On the one hand, that reduces the time available for the cloud to 
stretch the field lines (§ |4.6| ), therefore stopping the growth of the magnetic energy. On the other 
hand, since the cloud approaches more quickly the computational boundary, outflow of magnetic 
energy becomes soon important and significantly affects the growth curve of the magnetic energy. 

Recognizing these limitations we can still calculate values of the slope (m) and the intercept 
(A) for all cases with 9 > 30° , restricting ourselves to the growing portions of the curves for set 
1 and 3. The last four columns of Table |l] report values of the slope as well as of the intercept 
for the curves in Figure || of the magnetic energy in the total grid (m, ^4) and in the "cloud nose 
region" (m n , A n ) respectively, as obtained with a simple least chi — square method ( [Taylor 1982| ). 
For cases of set 1 and 3 the resulting slopes are somewhat sensitive to the choice of the "growing 
portion" of the curves. We point out that we report these values with the only purpose of giving 
a qualitative description of the evolution of the quantities considered; therefore they should not 
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be taken as to indicate anything too precise. Values of m and m n given in table average around 
in = 1.43 and in n = 0.9 respectively. In the following we shall simply take the representative 
values, m ~ 1.5 and m n ~ 1.. Finally, it is worth noting that for 9 > 30° the distance between the 
different curves (given by the difference in their A coefficients) is consistent with what expected 
from eq. 4-3 for different values of M and 9. This further supports the approach of eq. |4-3|. 



4.5. Dimensionality Issues 

In order to validate our 2D MHD numerical calculations, it is important to address the limits 
imposed by the symmetry along the cylinder axis due to the 2D geometry. For that purpose we 
carried out an analogous 3D version of Case 8, simulating the supersonic motion (M = 1.5) of 
a dense (x = 100), finite cylinder through a transverse magnetic field {9 = 90°). The length 
of the cylinder semi-axis was i z = 3R C , where R c was its radius. Motion was perpendicular to 
the cloud axis. Exploiting the mirror symmetry across the planes z = and y = 0, we were 
able to limit the calculation to only one quarter of the total actual volume. Therefore we used 
a computational domain x=[0,20], y=[0,8] and z=[0,8] with the cloud center initially located at 
x c = 2R C , y c = z c = 0. The resolution was 16 zones across a cloud radius with R c = 1 in numerical 
units. In order to make a sensible comparison with Case 8 we calculated the growth of the 
magnetic energy in the plane z = 0. In fact, if the 2D approximation is realistic, then the evolution 
of the magnetic energy on the plane z = resulting from the 3D calculation should resemble that 
seen in Case 8. The results from the 3D calculation are shown in Figure |(| In the three panels we 
report: log plots of the relative increment of the magnetic energy on the whole grid (top) and in 
the usual region of size R c x 2R C before the cloud nose (middle) versus time expressed in units 
of the sound crossing time . Also we report the evolution of the relative change of the kinetic 
energy (bottom). The Figure shows that up to about one r cr , the growth of the magnetic energy 



is qualitatively the same as in the 2D calculations. Following the same procedure as in §4.4, for 
this case we infer (m,A) = (1.26, 1.3) and (M n ,A n ) = (0.8,0.6). This slope is slightly smaller than 
in the previous case, meaning that in 3D, the growth of the magnetic field is somewhat slower, as 
expected, but is still faster than linear with time. 



4.6. Timescale 

Given the slope 1 ^ m ^ 1.5, in this section we derive the timescale T ma , to transform half of 
the cloud kinetic energy into magnetic form (formation of the magnetic shield). The cloud initial 
kinetic energy, Eko, can be expressed in terms of Ebo, the initial magnetic energy in the region of 
volume V s where field line stretching occurs, as 

Ew = \lxM 2 f3 ^ E B0 (4-5) 
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where V c is the volume of the cloud. Thus, setting the magnetic energy at t = T ma equal to half 



the initial cloud kinetic energy, and considering eq. 4-3 we obtain 



■tIXM 0— - 1 ~ -jxM (3— = — - ~ 2 M ; (4-6) 



4 ,/v V. 7 4 ,/v V, V r r 7 m\ 



(the numerical factor in the last expression derives from the assumption m ~ 1). From eq. 4-6 
derive 



we 



fiM] 



1 , , 2 -, . 1 

— S/f — — 1 



n M™ _i (V c 



V 8 / sinfl VVs 



r cs (4-7) 



and for m = 1.5 



using a representative sin# ~ sin 45° = l/y/2 and neglecting the factor (V c /V s ) 2 / 3 which probably 
amounts to a few tenths. This is the timescale over which a cloud gives up half of its kinetic 
energy to the background magnetic field by stretching the field lines. When viewed the other way 
around, r ma is also the timescale over which a cloud embedded in a magnetized flow is accelerated 
to move with the background plasma through magnetic forces. Therefore, r ma can be considered a 
"magnetic acceleration" timescale. We emphasize again such description is valid only up to several 



4.7. Physical Issues 

We CcLn now compare T ma with, the drag timescale T^ r (eq. |2-2[) , over which a cloud, is stopped, 
by the ram pressure (§@). Their ratio (using 7 = 5/3, 9 = 45° and /3 ~ a few) is: 

^ ~ X ^- X M 2 l m . (4-9) 

1~dr 



Since m > 1, Eq. 4-9 shows that for small Mach numbers (M) and large density contrasts (x) the 
magnetic tension is more important than the drag force. For example, for a "standard" interstellar 
cloud characterized by x ~ 100 and M ~ 1.5 (if m = 1.5) r ma /rd r ~ 0.37, so that a significant 
fraction of the cloud kinetic energy can be transformed into magnetic form. We point out that we 
use the drag time here instead, for example, of the crushing time because as long as the cloud is 
moving the magnetic field will be enhanced by stretching even after the passage of the crushing 
shock. Table § summarizes the relevant timescales characterizing the dynamical eveolution of a 
typical cloud in motion. 

In Figure ^ we report cuts of the log of the magnetic energy density (B 2 /2, in numerical units) 
along the axis of motion of the cloud for the following three representative cases: Case 3 (solid 
line), Case 7 (dotted line) and Case 11 (dashed line). The time corresponding to each of these cuts 
is expressed in terms of various dynamical timescales in Tab. 0. Values in column 6 show that 
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the evolutionary time is still below several r cr so that instabilities in the third dimension have not 



become too disruptive yet. In the same table we also calculate the ratio in eq. 4-9 for each case. 



The three lines in Figure |?j clearly exhibit a peak corresponding to the magnetic shield (§|2|). 
For Case 3, the magnetic energy density B 2 /2 ~ 40 is comparable with the ram pressure of the 
background flow (~ piv 2 ~ 100), in accord with Jones et al. (1996). For this case T ma /Td r ~9>1. 
Therefore, the cloud is mostly being decelerated by the ram presure of the flow and soon after, 
as the magnetic field lines relax, the magnetic energy drops. This effect, along with magnetic 
energy outflow through the boundary, is responsible for the turn-over visible in Figure [3| for Cases 
2 and 3. Similarly, for Case 11 (dashed line), since r ma ^> T^ r , only a small fraction of the cloud 
kinetic energy is transformed into magnetic energy. However, because of the large mass of the 
cloud (x S> 1), in this case even a small fraction of the cloud kinetic energy is enough to amplify 
the magnetic field in the magnetic shield beyond the background ram pressure limit. Finally, in 
Case 7 (dotted line), r ma /rrf r ~ 0.5 and, as we can see, B 2 /2 (~ 12.5) is much larger than the ram 
pressure limit (p ra , m = pv 2 ~ 2.25). The magnetic energy is still a small fraction (~ 20 %, see Fig. 
||) of the initial cloud kinetic energy, however. This is roughly in accord with the model []. 

We shall now return to the issue raised in §|2| about the limit on the magnetic energy that can 
be generated in the magnetic shield. It is clear that as long as the lines are being stretched by the 
gas motion, then the magnetic field energy will increase. This result was tested with a numerical 
simulation of a very large density (x = 10 3 ) 2D cloud with jet-like shape (in fact a sheet of dense 
gas). This object had a very large kinetic energy in its motion with respect to the background 
gas and it turned out that the magnetic energy grows (almost linearly) until it becomes so large 
on the cloud nose that the code cannot handle the calculation properly anymore. However, by 
contrast, it is usually observed in the solar system that there is a sort of "equilibrium" between 
ram pressure and magnetic pressure around massive objects moving supersonically with respect 
to the solar wind. This is the case for example at the earth magnetopause (e.g., Parks 1991| ) 



and also found in simulated interations between the solar wind and outflows from comets (Yi et 
al. 1996). That different result is due to the fact that in such cases the physical system reaches a 
stationary state, whereas the clouds considered here are never in a stationary flow. It is easy then 
to show that in a steady flow the previous limiting relation between magnetic pressure and ram 
pressure must hold. It is important to observe that stationarity is achieved in the aforementioned 
solar system cases because the basic motions of those objects are not affected by the flow and 
because the structural changes brought about through their interactions with the solar wind 
are strongly resisted. In fact the magnetopause is supported by a dipole magnetic field which 
responds strongly to any attempt of deformation (B 2 oc i?~ 6 ). Also the comet outflow provides a 
momentum flux (dynamic pressure) that varies as the inverse square of the distance (oc R~ 2 ). So 
again it is difficult to deform the shape of this kind of object. These facts allow the magnetic field 



Acording to eq. 4-3 (with m — 1.5), after a time t/r ma ~ 0.42 only a fraction ~ ji x (0.42) 1 ' 5 ~ 0.14 of the 



total initial kinetic energy is converted into magnetic form. 
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lines to easily slip around the surface for a "quick" attainment of stationary state. In our 2D cloud 
simulations the field lines are confined in the plane, so that they are bound to be stretched as long 
as the cloud moves. Growth is ultimately limited by deceleration of the clouds. However, even 
in 3D calculations since clouds undergo strong deformations on their leading edge the stationary 
state is never reached and magnetic flux can become trapped. In fact both the magnetic field and 
the cloud shape change continously in response to each other dynamical effects. Even in 3D MHD 
cloud simulations it turns out that the magnetic pressure in the shield becomes larger than the 
ram pressure by a factor of at least a few ( |Gregori et al. 1998| ). 

With this perspective, it can be concluded then that an accurate estimate of the upper limit 
for the magnetic energy in the magnetic shield can only come from appropriate 3D numerical 
simulations. However, it is immediately clear from the investigation presented in this paper that 
previous studies have underestimated the extent to which the magnetic field can be enhanced by 
stream motions in the ISM and overlooked the dynamical importance that this process can have 
in terms of the evolution of these streams as well as the energetics of the global ISM. In addition, 
as it will be properly addressed in the following section, this amplification of the magnetic field 
originating from the cloud motion has interesting observational implications. 

5. Comparison with Observations 

We attempt now to link our results to some recently found, interesting features of the 
magnetic field in diffuse clouds. As pointed out in §|3|, the conditions in the ISM are such that 
ideal MHD is applicable even to H I clouds characterized by the scales of interest here. 

Myers et al. (1995) compared data of 21 cm Zeeman effect obtained with the 26 m radio 
telescope at Hat Creek Radio Observatory (HCRO) and those with the 100 m Effelsberg telescope. 
They found, even though the magnetic strengths obtained with the two telescopes are usually 
consistent within experimental errors, that there are cases where the field strengths are apparently 
discrepant. For instance, in the particular direction (I, b)= (141.1, 38.8), the HCRO gives 
Bn = 18.9 ± 1.8 fj,G, whereas the 100 m measurement is 3.5 ± 3.7 [iG. Those authors argued 
that this is probably real and due to the significant variation of the magnetic field strength inside 
the larger HCRO beam compared to that in the 100 m beam. This means that, in that region, the 
field has a significant structure on a length scale L < 1.6pc and measuraments of B\\ > 18.9/xG 
should be found by the 100 m telescope somewhere within the HCRO beam. In §|2| we have 
seen that a cloud propagation through a magnetized background is able to modify significantly 
the initially uniform field structure. This is particularly true for 9 > 30°. As already seen in 
Figure 0, after a time of the order of r cr , the intensity of the field in the magnetic shield can be 
several times as large as the background value. The magnetic shield dimensions are of the order 
of the cloud radius (R c ). The variation in the magnetic field measured by Myers et al. (1995) 
(from 3.5 ± 3.7 ^G up to values > 18.9) is consistent with the amplification resulting from our 
numerical study. Therefore, it is plausible to explain the local inhomogeneities observed in Myers 
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et al. (1995) in terms of the interaction between the cloud motion and the background magnetic 
field. 

Another interesting observational finding is that in many H I clouds, the magnetic and kinetic 
energy densities have comparable values (Heiles 1989, Vershuur 1989, Myers & Khersonsky 1995, 
Myers et al. 1995). For some cases the magnetic energy is actually in excess of the kinetic (Myers 
et al. 1995). Myers et al. (1995) point out that if collisions between neutral and ions are high 
enough, MHD waves will propagate through the field-fluid system at the Alfven speed based on 
the combined neutral/ionic mass density. The kinetic and magnetic energy densities in an Alfven 
wave are in equipartition. Therefore, if the magnetic energy in the waves is similar to that in the 
mean magnetic field, the nonthermal line broadening would correspond to the Alfven speed for the 
neutral/ionic mass. In other words, the measured magnetic and kinetic energies of the gas should 
be similar. On the other hand, those observational results are in agreement with the depiction 
given in this paper, according to which the magnetic energy is enhanced by field lines stretching 
produced by coherent gas flow, i.e., by cloud motions. Since this field enhancement occurs at the 
expense of the cloud kinetic energy, it makes sense that some equipartition between magnetic and 
kinetic energy is achieved. Eventually the cloud is stopped, either by the magnetic field tension 
or by the ram pressure of the background flow, so that at some point the kinetic energy becomes 
even smaller than the magnetic energy (although the latter must still be less than the initial cloud 
kinetic energy). On the other hand we have seen that only a weak amplification of the magnetic 
field occurs, when the latter is mainly aligned with the cloud motion. Therefore, not for all clouds 
magnetic and kinetic energy are expected to be comparable. 

6. Conclusions & Summary 

We have studied the motion of an individual 2D clouds through a magnetized lower density 
background medium, for different density contrasts, x> Mach numbers, M, and field geometries, 9. 
Our major findings can be summarized in the following items: 

• A cloud moving through a magnetized medium tends to interact with the magnetic field 
in a way that depends upon the angle 6 between the initial cloud velocity and background 
magnetic field. According to our simulations, for large angles 6 the magnetic field lines are 
stretched efficiently by the cloud motion and, therefore, the magnetic field is significantly 
amplified at the expense of the cloud kinetic energy. This is demonstrated here for angles 
6 > 30°. Our simulations then suggest that when 9 > 30° the magnetic shield developed by 
a moving cloud may behave as an efficient "bumper" in MHD cloud collisions ( |Miniati~e^ 
al. 199fl ). 

• Our numerical experiments show a global magnetic energy growth reasonably well described 
by a power law in time with index m ~ 1.2—1.5 for a wide range of values of density 
contrast (x), Mach number (M) and field geometry (9). This behavior is not significantly 
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affected by the presence of radiative losses. In addition, it seems confirmed by at least 
one analogous 3D simulation as well (m ~ 1.2,1.3). We define a characteristic timescale 

J_ 2-m 

T ma ~ \PX) m M rn 4^ for part of the cloud kinetic energy to be transformed into magnetic 
energy. We note, as well, that r ma can also be interpreted as the time for a cloud to be 
accelerated to the flow velocity by Maxwell stresses. Since m > 1, for high x, low Mach 
number clouds this dynamical timescale would be shorter and therefore more relevant then 

m— 1 2_ 

the drag timescale r^ r ~ x m M m T ma . 

• We suggest that this mechanism, if confirmed more generally by 3D MHD simulations, 
could provide an explanation for some characteristics observed in the interstellar magnetic 
field. In particular the existence of magnetic field inhomogeneities on a scale about 1 pc in 
correspondence of H I clouds can be interpreted as the formation of the magnetic shield, as 
outlined in §|B[ In addition the equipartition of kinetic and magnetic energy, independent of 
whether or not the cloud is self-gravitating, would follow naturally from this depiction where 
the the magnetic energy is provided by the cloud kinetic energy itself. 

• In the majority of the cases we simulated the magnetic field lines drape around the 
propagating cloud, even though no strong magnetic shield forms. In addition, as the 
symmetry across the motion axis is broken (for 0° < 8 < 90°), the flow in the wake becomes 
more turbulent, thus facilitating the onset of tearing mode instabilities which weaken the 
magnetic field strength there. This means that the "flux rope" mentioned by earlier authors 
might develop only in very special cases or exist only temporarily. Asymmetries in the 
magnetic tension around the cloud can induce vortical motions within the cloud itself. Those 
motions, in turn, could influence its subsequent development or annihilation of fields within 
the cloud. 

We are grateful to Ruth Dgani, the referee, for insightful comments and interesting discussions 
which improved the manuscript. We acknowledge the support provided by Gianluca Gregori in 
performing the 3D MHD numerical calculation presented in this paper. The work by FM and 
TWJ was supported in part by the NSF through grants AST-9619438, INT-9511654, by Nasa 
grant NAG5-5055 and by the University of Minnesota Supercomputer Institute. The work by DR 
was supported in part by KOSEF through grants 975-0200-006-2 and 981-0203-011-2. 

REFERENCES 

Chandrasekhar, S. 1961, Hydrodynamic and Hydromagnetic Stability (Oxford: Claredon) 
Dgani, R. 1998, Rev. Mexicana Astron. Astrophis., 7, 149 
Dgani, R., & Soker, N. 1998, ApJ, 495, 337 



-19- 



Draine, B. T., k McKee, C. F. 1993, ARA&A, 31, 373 

Ferrara, A., k Field, G. B. 1994, ApJ, 423, 665 

Gardner, F., k Whiteoak, J. B. 1963, Nature, 197, 1162 

Gregori, G., Miniati, F., Ryu, D., k Jones, T. W. in preparation 

Gomez De Castro, A. I., Pudritz, R. E., k Bastien, P. 1997, ApJ476, 717 

Heiles, C. 1976, ARA&A, 14, 1 

Heiles, C. 1989, ApJ, 336, 808 

Klein, R. L, McKee, C. F., & Colella, P. 1994, ApJ, 420, 213 
Jones, T. W., Kang, H., k Tregillis I. L. 1994, ApJ, 432, 194 
Jones, T. W., Ryu, D., k Tregillis I. L. 1996, ApJ, 473, 365 

Mac Low, M., McKee, C.F., Klein, R.I., Stone, J.M., k Norman, M.L. 1994, ApJ, 433, 757 

Malagoli, A., Bodo, G. k Rosner, S. 1996, ApJ, 456, 708 

Mathewson, D. S., k Ford, V. L. 1970, MNRAS, 74, 139 

Miniati, F., Jones, T. W., Ferrara, A., k Ryu, D. 1997, ApJ, 491, 216 

Miniati, F., Ryu, D., Ferrara, A., k Jones, T. W. 1998, ApJ, in Press 

Morris, D., k Berge, G. 1964, ApJ, 139, 1388 

Murray, S.D., White, S.D.M., Blondin, J.M., k Lin, D.N.C. 1993, ApJ, 407, 588 
Myers, P. C., Goodman, A. A., Glisten, R., k Heiles, C. 1995, ApJ, 442, 177 
Myers, P. C., k Khersonsky, V. K. 1995, ApJ, 442, 186 

Parks, G. K. 1991, Physics of Space Plasmas: An Introduction, (Addison Wesley) 

Ryu, D., k Jones, T. W. 1995, ApJ, 442, 228 

Ryu, D., Jones, T. W. k Frank, A. 1995, ApJ, 452, 785 

Ryu, D., Miniati, F., Jones, T. W., k Frank, A., 1998, ApJ, in Press 

Schiano, A. V. R., Christiansen, W. A. k Knerr, J. M. 1995, ApJ, 439, 237 

Soker, N., k Dgani, R. 1997, ApJ, 484, 277 

Spitzer, L. 1978, Physical Processes in the Interstellar Medium, (New York: Wiley-Interscience) 



-20- 



Spoelstra, T. A. Th. 1977, Usp. Fiz. (SSSR) 121, 336 

Taylor, J. R. 1982, An Introduction to Error Analysis, The Study of Uncertainties in Physical 
Measurements, (University Science Books) 

Valee, J. P., & Kronberg, P. P. 1975, A&A, 43 233 

Verschuur, G. L. 1970, "Observational Aspects of the Galactic Magnetic Field" in Interstellar Gas 
Dynamics (Dordrecht: D. Reidel) 

Verschuur, G. L. 1989, ApJ, 339, 163 

Vietri, M., Ferrara, A., & Miniati, F. 1997, ApJ, 483, 262 

Yi, Y., Walker, R. J., Ogino, T., Brandt, J. C. 1996, J. Geophys. Res., 101, 27585 

Zeldovich, Ya.B., Ruzmaikin, A.A.,& Sokoloff, D.D. 1983, Magnetic Fields in Astrophysics, (New 
York: Gordon and Breach) 



This preprint was prepared with the AAS IATeX macros v4.0. 



-21 - 



Table 1. Summary of 2D-MHD Single Cloud Simulations 











0= 4, 


7=5/3 










CASE a 


X 


M b 




(°) 


Grid-Size c 
(^) 


End Time d 

(Tcs) 


m e 


A e 




A f 

Sin 



l 


10 


10 





10 x 5 


4.2 










2 


10 


10 


45 


20 x 10 


6.6 


1.4 


0.2 


0.8 


1.6 


3 


10 


10 


90 


10 x 5 


4.2 


1.5 


0.7 


0.6 


1.8 


4 


100 


1.5 





15 x 15 


13 










5 


100 


1.5 


10 


20 x 10 


15 










6 


100 


1.5 


30 


20 x 10 


15 


1.4 


-1.4 


1.1 


0.1 


7 


100 


1.5 


45 


20 x 10 


13 


1.4 


-1.2 


1.0 


0.4 


8 


100 


1.5 


90 


20 x 7.5 


9 


1.5 


-1.3 


1.0 


0.6 


9 


100 


1.5 


90 


20 x 7.5 


9 


1.4 


-1.2 


0.9 


0.7 


10 


10 


1.5 


45 


20 x 10 


9.6 


1.2 


-1.3 


0.6 


0.3 


11 


100 


10 


45 


20 x 10 


9.6 


1.6 


0.3 


1.2 


1.8 



a All models use (3 = 4, 7 = 5/3. 

b The Mach number is referred to the intercloud sound speed, c s i- 

c The grid size is expressed in units of cloud radius. One cloud radius has i? c =50 zones. 
d The end time is expressed in terms of sound crossing times t cs . 

Coefficients for the growth of the magnetic energy over the whole grid. For cases 2, 3, 9 and 10 only the growing 
portion of the magnetic energy growth curves has been used to determine m and A (see the text for details). These 
are not provided for Cases 1, 4 and 5 because the model of § 4.4 is not appropriate for these cases. 

Coefficients for the growth of the magnetic energy inside a region of size R c x 2R C , placed before the cloud's nose. 
The same precautions as in the previous note apply here. 



Table 2. Timescales 



CASE 


X 


M 


Tc 1 Tcs 


Tcr } Tcs 


(Tma/Tcs) SU1 9* 




1,2,3 


10 


10 


0.1 


0.63 


8.8 


l 


4,5,6,7,8,9 


100 


1.5 


0.67 


13.3 


21.8 


66.7 


10 


10 


1.5 


0.67 


4.2 


4.7 


6.67 


11 


100 


10 


0.1 


2 


41 


10 



a The factor sin# is to take into account the slight difference between cases characterized by the same \ an d M, 
but different angles 9. For Cases 1, 4 and 5 (9 < 30°), as already mentioned in §1.4, r ma is meaningless. 
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Table 3. Cloud Timescales 



CASE 



3 8.8 0.3 2.4 24 3.8 

7 0.46 0.42 0.19 19.5 0.95 

11 5.8 0.14 0.84 84 4.2 
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FIGURE CAPTIONS 



Fig. 1. — Diagram of a cloud moving through a uniformly magnetized background medium with a 
velocity making an angle 9 with the field lines. 

Fig. 2. — Density distribution with superimposed magnetic field lines for a single cloud propagating 
through a magnetized medium. Field lines are contours of the magnetic flux. Density figures are 
inverted grayscale images of log(p). In both cases (3 = 4, x = 100 an d M = 1.5. Top panel, for 
9 = 45° and t = 13r cs = 19.5r c , exhibits the formation of a consistent magnetic shield. Bottom 
panel is 9 = 10° and t = 15t cs = 22.5r c ; the field lines are draping around the cloud, although the 
magnetic shield is not as developed as in the previous case. 

Fig. 3. — Time evolution of the magnetic energy enhancement (AEb/Exo, left panels) and of the 
kinetic energy (Ek/Ekq, right panels), normalized to the initial cloud kinetic energy, Eko, for 
individual clouds propagating through a magnetized medium. Time is espressed in units of the 
background sound crossing time (t cs = Rc/c s i) for all cases. Time labels, however, differ because 
the timescales for the evolution of magnetic and kinetic energy strongly depend on the parameters 
(x and M) defining the various cases (see §||| and §4^ for details). In the top panels (set 1) open 



triangles correspond to Case 1, open circles to Case 2 and filled circles to Case 3 respectively. In the 
middle panels (set 2) open squares represent Case 4, filled squares Case 5, open triangles Case 6, 
open circles Case 7, filled circles Case 8 and filled triangles Case 9. In the bottom panels open and 
filled circles correspond to Case 10 and 11 respectively. The middle, left panel shows that for Case 
4 (9 = 0°, open squares) the magnetic energy growth is slightly larger than for Case 5 (9 = 10°, 
filled squares). This is only due to the choice of the normalization parameter Eko which, due to 
the setting of the grid, is twice as large for Case 5 than for Case 4. Compare with Figure || for 
the relative growth of the magnetic energy normalized to its initial value. Finally, beware that the 
sudden drop in the kinetic energy for Case 9 (filled triangle) is a numerical rather than a physical 
effect. 

Fig. 4. — Magnetic field geometry for Case 5 (x = 100, M = 1.5,9 = 10°) at three different times 
corresponding, from top to bottom, to: t/r cs = t/(Mr c ) = 12, 13, 14. As in Fig. E| field lines are 
contours of the magnetic flux function. Finer contouring is used, however, with respect to Fig. [2|. 
The sequence of panels displays the evolution of the "tearing-mode" instability in the cloud wake, 
with the formation of closed field line loops through reconnection and their subsequent annihilation. 



Fig. 5. — Time evolution of the relative magnetic energy enhancement (AEb/Ebo) relative to the 
whole numerical grid (left panels) and to a region of size R c x 2R C placed before the cloud nose 
(right panels). In each case Ebo is the initial magnetic energy relative to the region considered. 
Curves correspond to log plots of AEb/Ebo versus time normalized to the sound crossing time t cs . 
As in Figure |3], time labels change for the various cases, according to the different characteristic 
timescales. Each case of Table [l] is represented in the same panel and with the same type of point 
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as in Figure |3|. Note that for a better representation of Case 2 and 3 the less interesting Case 1 
was not plotted in the top, left panel. 

Fig. 6. — 3D analogous of Case 8 in Table [l]. For this case \ = 100, M = 1.5 and 6 = 90°. Top and 
middle panels correspond to the filled circle curves in the middle, left and right panels of Figure [| 
respectively. Bottom panel, on the other hand, corresponds to the filled circle plot of the middle, 
right panel of Figure [3|. 

Fig. 7. — Cuts through the cloud and along the cloud axis of motion, of the log of the magnetic 
energy density (B 2 /2) expressed in numerical units. Solid line corresponds to Case 3, dotted line 
to Case 7 and dashed line to Case 10. Times corresponding to each curve and expressed in different 
units are reported in tab. ||. In all three cases it is evident a peak in the magnetic energy density 
due to the presence of the magnetic shield in front of the cloud. 



This figure "fig2.jpg" is available in "jpg" format from: 
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This figure "fig4.jpg" is available in "jpg" format from: 
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